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We present ah initio calculations for single-atom vacancies in graphene. In agreement with earlier 
work, we find that vacancies are responsible for the magnetism recently observed experimentally, 
creating a center with net spin. For small supercells, there is a strong symmetry breaking in the 
supercell with respect to the perfect graphene structure, and this symmetry breaking is much smaller 
for larger supercells. The influence of adsorption of H, O and N atoms on the spin center of the bare 
single- atom vacancy and its influence on the magnetism induced in these samples is investigated. A 
rich variety of mechanisms is found for the generation of magnetic moments in this system. 

PACS numbers: 81.05.Uw, 61.72.Ji, 73.22.-f 



Macroscopic magnetic ordering phenomena in organic 
materials has been one of the most exciting and interest- 
ing subjects in physics, chemistry and materials science. 
Organic magnetic materials have attracted the interest of 
scientists not only because of their fundamental proper- 
ties Q , but mainly because of their potential applications 
in high-tech devices. Among them, we have to mention 
applications in medicine and biology, where organic mag- 
netic materials would be highly desirable due to inherent 
biocompatibility 3]. 

Magnetism has been observed recently in graphite sam- 
ples ^ H, 1^ and its origin is attributed to the creation 
of vacancies in the material, which would lead to the 
appearance of unpaired spins and consequent magnetic 
moments. Previous studies (Jj, Q have confirmed that 
single-atom vacancies indeed lead to the appearance of 
a center with a magnetic moment. A recent study has 
shown that the existence of magnetism in vacancies in- 
duced by the removal of more than one atom can lead to 
non-magnetic states, as well as magnetic ones 0. 

Vacancy creation in graphite samples can be achieved 
either by ion irradiation 7] or chemical modification 
After almost two years of the production of magnetic 
graphite samples by some of the present authors using 
the chemical modification method its magnetism re- 
mains just as strong as before, at room temperature and 
air exposure. It is natural, then, to ask what happens 
when a foreign species is adsorbed on the magnetic cen- 
ter. In this Letter, we analyze how the magnetic moment 
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induced by single-atom vacancies in graphene (SAV-g) 
changes when a single atom of H, N or O is adsorbed on 
the magnetic center. As will be shown, the adsorption 
of these atomic species changes the fundamental mecha- 
nisms through which the magnetic moment is created in 
very different and surprising ways. 

Our calculations have been performed in the frame- 
work of Densit y Fu nctional Theory [Tol . ITU using the pro- 
gram SIESTA hj. We have used Troullier-Martins pseu- 
dopotentials 13], including nonlinear core corrections for 
N, the PBE-GGA exchange-correlation functional 
and a split- valence double-zeta basis set with polarization 
orbitals (DZP) for all atoms. Because the Fermi surface 
of graphene is a single point, Monkhorst-Pack grids 
of 9x9x1 and 5x5x1 points were used for atomic relax- 
ations in the small (4x4x1, SSC) and large (8x8x1, LSC) 
supercells, respectively, which amount to equivalent k- 
point samplings of 36x36x1 and 40x40x1 in the Brilluoin 
Zone (BZ) of a unit cell of graphene. Densities of states 
were calculated with denser grids, equivalent to BZ sam- 
plings of 96x96x3 and a gaussian smearing of 0.08 eV. 
The height of the supercells was fixed to c/a = 4.2 in 
the case of SSC, and c/a = 8.2 for LSC. All forces on 
atoms were smaller than 0.02 eV/A, and cell relaxations 
included an additional constraint requiring that pressures 
be smaller than 0.05 GPa. Throughout the text, all net 
spin charges, = (j| — cr|, are referred to in units of 
the electron spin, that is, (jg = |, and net spin densities, 
Spcr{r) = p]{r) — Pii^) are in units of (5<j/Bohr^. 

For both SSCs and LSCs, simple relaxation of atomic 
coordinates with of bare SAV-g, simple relaxation of 
atomic coordinates with supercell vectors parallel to 
those of a unit cell of perfect graphene yields stress ten- 
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sors with non-zero xy and yx elements with magnitude 
comparable to that of the xx component; furthermore, 
the XX and yy components have very different magni- 
tudes. Therefore, there is room for further relaxation of 
the supercell. Variable-cell relaxation yields lattice vec- 
tors with |ai| 7^ \d2\ and = 61.5° for a bare SAV-g 
SSC. For the bare SAV-g LSC, we have |ai| = \a2\ and 
= 59.5°, which shows that the in-plane interaction be- 
tween vacancies plays an important role in the process of 
changes of supercell symmetry. The spin-dependent total 
density of states (TDOS) per C atom has two peaks close 
to the Fermi level, as shown in Figs, ^a) and (b), but 
it is clear that the peaks for the SSC TDOS are higher 
and farther from each other than the ones in the LSC 
TDOS. It can also be seen from the 2p projected den- 
sity of states (PDOS) per C atom (Fig. H^c)) that the 
TDOS features are predominantly determined by those 
states. Therefore, in agreement with earlier DFT calcu- 
lations ^9] , we find 2p states to be the responsible by the 
observed magnetism in bare SAV-g. 
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FIG. 1: (Color online) (a) Spin-dependent TDOS, per atom, 
for SSC of bare SAV-g. (b) Spin-dependent PDOS of all C 
atoms of SSC with bare SAV-g. (c) TDOS of LSC of bare 
SAV-g. Red, solid lines: spin up. Blue, dashed lines: spin 
down. 



Denoting the SAV-g with a certain adsorbed atomic 
species X by X-SAV-g, we start by analyzing H adsorp- 
tion on the magnetic center of a SAV-g cell. Both the 
SSC and the LSC are allowed for complete relaxation, 
undergoing a similar symmetry breaking of the unit cell 
as in the case for SAV-g. The net magnetic moment is 
zero, in agreement with the results of Lehtinen et al 7]. 
The H atom undergoes an out-of-plane displacement, and 
for the SSC the overall distortion of the sheet is very 
pronounced, while it is smaller in the LSC. This sug- 
gests that relief of mechanical tensions by cell-symmetry 
breaking are rather important, which and for the SSC of 
a H-SAV-g with for supercell vectors parallel to the ones 
of the perfect graphene unit cell, the stress tensors also 
have nonzero xy and yx components, as well as different 
XX and yy elements. However, the surprising fact is that 



under applied mechanical tension, a liquid spin appears 
in both supercells, being 6a =0.94 and 0.25 for the SSC 
and the LSC, respectively. This fact suggests that me- 
chanical tension is indeed the main responsible for this 
different spin, because in the LSC, there is more room for 
tension distribution. Fig. [21 shows the relaxed structure 
and the net spin charge density for the H-SAV-g LSC. 
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FIG. 2: (Color online) (a) Relaxed SSC of H-SAV-g. (b) 
Liquid spin density of SSC for supercell vectors parallel to 
those of a perfect graphene sheet. 

While a more detailed determination of this mecha- 
nism demands a deeper investigation, an analysis of the 
PDOS and spin-difference charge density distribution, 
5pcr{r)^ shows that it involves the polarization of C atoms, 
while the H states contribute with no spin polarization. 
The spin-dependent 2p PDOS per C atom, as depicted 
in Fig. OJb), clearly shows that the spin polarization is 
completely onto the atoms of this species, while the Is 
states of the H atom, shown in Fig. OIJc), do not con- 
tribute to the magnetism in a significant way, since its 
contribution is about four times smaller than the aver- 
age contribution of the C atoms. Moreover, in contrast 
to the case of bare SAV-g, in which the magnetic mo- 
ment is concentrated mostly on the C atom onto which 
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the foreign atom is adsorbed, large magnetic moments 
also on the two other atoms which approach to form a 
pentagon. Even more interesting is the fact that the or- 
dering of the spin polarization on the surrounding atoms 
has an alternating disposition, but with decaying inten- 
sity as we move further away from the vacancy. Since 
the system exhibits a net magnetic moment due to the 
imbalance between spin up and spin down polarizations 
this arrangement is ferrimagnetic. 



vacancies, the N atom will not enter substitutionally and 
such a configuration is more likely to occur. 
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FIG. 3: (Color online) (a) Total spin-up and down DOS, per 
atom of LSC of H-SAV-g, shown only in the vicinity of the 
Fermi energy (set to zero), where the two DOS differ more 
significantly from each other, (b) Spin-up and spin-down 2p 
PDOS, per C atom, in a LSC of H-sav-g. (c) PDOS of IS 
states of the adsorbed H atom in the LSC of H-sav-g. Red, 
solid lines: spin up. Blue, dashed lines: spin down. 



Next, we analyze the influence of N adsorption onto 
the magnetic center of a bare SAV-g, in the geometry 
depicted in Fig. E|^a). Both SSC and LSC, whether 
cell-relaxed or not, display Spa- = 1, that is, the mag- 
netic moment in N-SAV-g is independent of geometri- 
cal details. To better understand this, we plot spin- 
difference charge density isosurfaces with different values 
to check how liquid spin is distrbuted in the supercell. 
For a spin-difference of Spa- = ±0.01, we can see that 
the positive spin charge is widely distributed around the 
adsorbed N atom, while the negative spin charge is nar- 
rowly distributed around the closest C atom. Isosurfaces 
of Spa- = 0.2 are tightly bound to the N atom, while 
on C atoms Spa- = 0. Looking at the C-N bond length, 
ac-N = 1-25 A, we can infer that the N atom makes a 
double bond with the C atom. This causes the satura- 
tion of the unpaired C spin, with 5(j = 1 coming from 
the N atom. Since the net spin is tightly bound to the 
N atom, this explains the insensitivity of the magnetism 
to geometrical details, in contrast with the H-SAV-g. Al- 
though according to our calculations the ground state 
for a N atom would be that in which it enters the sheet 
substitutionally, the situation depicted in Fig. El could 
be found, experimentally, in a zigzag state for a vacancy 
where many atoms have been removed. In these large 



FIG. 4: (Color online) (a) Isosurface of 6pa = ±0.01 in SSC 
with N-SAV-g; positive 6 pa are depicted in red, negative ones 
in blue. The widest surface is around the N atom, and C 
atoms display a far smaller isosurface. (b) Isosurface of ^pcr = 
0.2. No negative values are observed on atoms other than C. 



The case of 0-adsorption is somewhat the opposite 
of N-adsorption. 0-adsorption on the magnetic center 
causes the magnetic moment to completely vanish, no 
matter if the cell is distorted or not, for both SSC and 
LSC. Although this is in striking constrast with the re- 
sults of Esquinazi et al jQ] , where the authors find that 
the presence of O2 triggers a magnetic moment, it is im- 
portant to bear in mind that the present calculations do 
not correspond to the experimental situation, in which 
the authors use molecular oxygen. Further theoretical 
work, with systems which approximate more realistically 
the experimental situation - that is, layered graphite with 
adsorbed and perhaps interstitial O2 as well - is neces- 
sary to clarify the situation. Fig. [51 shows the TDOS, 
per atom of the LSC, of 0-SAV-g. 

Our findings certainly will have important implica- 
tions in most of the applications of carbon-based mag- 
netic materials in nanotechnology, such as sensors, detec- 
tors, actuators, telecommunications, electronics, biosen- 
sors, magnetic materials separation, etc.. However, those 
implications will be particularly important for applica- 
tions in medicine, specially for targeted drug delivery by 
use of magnetic field gradients llfl • Our calculations 
imply that, if magnetic graphene is to be used in this kind 
of application, one must avoid the saturation of the mag- 
netic centers, as much as possible, by oxygen, and that a 
good situation might be that in which the magnetic cen- 
ters have N adsorbed. Experimentally, it is known that 
biologically relevant molecules such as cholesterol oxidase 
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FIG. 5: (Color online) Total density of states per atom of 
0-SAV-g LSC. 



with vacancies on different layers, and in different con- 
centrations, also exhibit a net magnetic moment 17]. 

In summary, we have studied three very different mech- 
anisms that lead to the appearance of a net magnetic 
moment on a graphene sheet with a single-atom vacancy, 
both bare and with foreign atomic species adsorbed on 
top of it, particularly H, N and O. Both small and large 
supercells have been studied. A strong breaking of the 
super cell symmetry is observed for all systems. While 
the magnetic moment initially present in the bare SAV- 
g, regardless of the geometrical details of the supercell, is 
completely quenched by O adsorption and is remarkably 
stable upon N adsorption, we find that it can be induced 
by tensile stress of the lattice when H is adsorbed in the 
system. Although we have studied only the adsorption 
of atomic species on the system, it is highly possible that 
adsorption of molecules is also capable of preserving the 
magnetism in the system, and its implications for device 
manufacturing make this a good subject for further stud- 
ies. 



attach to magnetic Fe304 nanoparticles through amino 
and carbodiimide-activated carboxyl groups 16], both of 
which contain N. This suggests that proteins could attach 
to magnetic graphene without the use of a ligand, which 
would make the use of magnetic graphene extremely at- 
tractive for biomedical applications. Our results also im- 
ply that, if graphene with vacancies should be used for 
any kind of application involving its magnetic proper- 
ties, it should be, as much as possible, oxygen- free. If 
H should be present in the system, then deposition of 
the graphene sheet onto a substrate might provide the 
tension necessary for the appearance of a magnetic mo- 
ment due to lattice mismatch. The use of multi-layered 
graphitic structures is, in principle, also not discarded: 
preliminary FPLAPW results show that bulk graphite 
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